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Abstract: 

We calculate the spin-averaged J/^-nucleon scattering length aj/^p by directly applying 
QCD sum rule to J /ip-N forward scattering amplitude. Our result aj^ = —0.10 ± 0.02 fm, 
predicts the possibility of bound states with nuclei, though the force is weaker than that of 
the light vector mesons (p, u>, <f>)-N cases. Up to dimension-4 gluonic operators, we evaluate 
the scattering length with twist- 2 contribution, which increases the absolute value of the 
scattering length about 30%. If we apply aj/^p to the effective mass of J ftp in nuclear matter 
on the basis of the linear density approximation, it shows very slight decrease (4 ~ 7 MeV) 
at normal matter density. 
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1 Introduction 



Theoretical analysis on the in-medium properties of hadrons is increasingly required by 
various on-going and forthcoming heavy-ion experiments (such as SPS, LHC (CERN) and 
AGS, RHIC (BNL)) [|TJ. In particular, experimentally it is important to observe vector 
mesons, because they decay into lepton pairs and carry the information inside the mat- 
ter without disturbance of the strong interaction. The properties of light vector mesons 
in nuclear matter have been extensively studied in various theoretical approaches such as 
effective hadronic models @ and QCD sum rules (QSR's) ||, f|, |5| f|. Vacuum properties of 
the vector mesons have been successfully studied by the QSR's |7|, ||]. The method enables 
us to express physical quantities such as mass and decay width in terms of the parameters 
in QCD Lagrangian and vacuum condensates. Extending the vacuum QSR to finite density, 
we can consistently incorporate the effects of nuclear matter into the form of in-medium 
condensates. There are two methodologically different ways for in-medium QSR. Firstly, 
T. Hatsuda and S.H. Lee developed the in- medium QSR formalism for light vector mesons 
||. They found a 10 ~ 20 % decrease of the p and u mesons at normal matter density. 
Secondly, for light vector mesons we formulated in-medium QSR j|, |6| based on the relation 
between a scattering length and a mass shift ||. In this approach with the Fermi gas model, 
in-medium correlation function is devided into vacuum part and one nucleon part. This one 
nucleon part corresponds to the forward vector meson-nucleon scattering amplitude. The 
QSR analysis on the forward scattering amplitude enables us to obtain the information 
for vector meson-nucleon interaction. Moreover, from the information we can estimate the 
change of spectra for vector mesons in nuclear matter. The difference between these two 
approaches has been discussed in |7j, ||. Eventually we derived in || that both of them are 
based on almost the same idea and can lead the consistent results with those of the effective 
models. 

In this paper we apply the QSR analysis established in || to a heavy quark system 
with equal mass for quark and antiquark. As a concrete system we focus on J/if>, which 
is a low-lying charmonium state ( 3 Si). To study the medium modification of J /if) has the 
following reasons: 

1. We have the detailed experimental information for the charmonium. In particular, 
the spectrum of J /if) is extremely narrow for the leptonic decay (IVj- ~ 5 keV). So it 
would be good tool to observe the change of the spectra (e.g., mass shift) in nuclear 
matter. 

2. Since charmonium and nucleons consist of quarks in different kinds of flavors, J/ip-N 
interaction is purely gluonic without quark exchange to first order in elastic scattering. 
This simplification reduces our practical calculation. 

3. Theoretical studies for J /if) in QSR have succeeded only in the description of the free 



state |R ITT 
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4. To utilize J /if) suppression [15] as a direct signal of the quark-gluon plasma (QGP) 



phase, we need to estimate the effect of nuclear absorption theoretically fL3f . It seems 



prompt to conclude that the present experimental data [Q can be explained only 
by nuclear absorption, until we investigate the J/^-nucleon interaction in detail. For 
that purpose, it is reasonable as the first step to study the J/ip-N elastic scattering 
at low energy. 

Motivated by these, we calculate J/ip-N scattering length and the mass shift of J /if} in 
nuclear matter. That is, the first aim is to estimate the essential features of the interaction 
between J/tf> and N through the scattering length. In practice, by applying QSR to J/ifi-N 
forward scattering amplitude we calculate the scattering length. The scattering length is 
a physically very important quantity in free space, because it is the unique observable in 
J/ifj-N elastic scattering at low energy. If it is negative, then we could predict attractive 
nuclear force capable of binding J/ if> to nucleus, so that J/ if> could lead to a bound state 
with nucleus. Moreover, due to absence of Pauli blocking unlike the case of light vector 
mesons- iV system, the effective J/if)-N interaction will not have a short-range repulsion. 
The prediction of such exotic state will give exciting new directions in nuclear physics. As 
is well known, since the nuclear force is repulsive for isovector mesons, n meson forms tc- 
nucleus bound state by Coulomb attractive force. On the other hand J /if) is expected to 
be bounded only by attractive nuclear force from the isoscalar property. As was pointed 
out in [|1|, this interaction should be sufficiently attractive to allow a bound state. The 



probability of such exotic states has recently discussed for rj, u and D cases [Tj| . 

The second aim is how the superposition of elementary J/ if)-N scattering at low energy 
affects the effective mass of J /if} in nuclear matter. When we work in the dilute nucleon 
gas, we find that the mass shift is linearly dependent on the density (linear density approx- 
imation). 

This paper is organized as follows. In section 2 we summarize the relation between the 
scattering length and the mass shift in the linear density approximation [Q. In the actual 
calculation we adopt moment sum rule method to the forward scattering amplitude. In 
section 3 the Wilson coefficients in the OPE side is explicitly given for twist- 2 operator. 
In section 4 in order to obtain unknown hadronic parameters for the forward scattering 
amplitude we apply moment sum rule to vacuum correlation function. In section 5 the 
numerical results of the scattering length and the mass shift of J /if) is shown. Finally 
concluding remarks are given. 



2 The relation between scattering length and mass 
shift 

Let us first review the relation between the scattering length and the mass shift on the basis 
of QSR method @, ||. The starting point of this approach is the following vector current 
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correlation function in the ground state of nuclear matter with nucleon density pn- 



K M (l) = i J d 4 xe^(TJ,(x)Jm) NMiPNh (2.1) 

where = (u, q) is the four-momentum carried by the J ftp vector meson current J M (ar) = 
c7 M c(a;) with the quantum numbers. Following the QSR method, when we apply an operator 
product expansion (OPE) to this correlator at deep Euclidean region (Q 2 = —q 2 > 0), it 
is supposed that the /^-dependence of this correlator is entirely contained into the p^- 
dependence of various condensates. Moreover we assume the Fermi gas model taking into 
account the Pauli principle among uncorrelated nucleons for the nuclear matter. In this 
approximation, in-medium correlation function reads 

/VF fj^Ti 
727v^ T - (g) ' (2 - 2) 
spin,isospiri \ I "0 

where 11° (q) is in- vacuum correlation function and J2 sp in,isos P in denotes the sum of spin and 
isospin states for nucleons in nuclear matter. T liv {q) is the vector current-nucleon forward 
scattering amplitude defined as 

T^(u,q) = i J d 4 xe iq - x (N( P s)\TJ,(x)4(Q)\N(ps)). (2.3) 

Here \N(ps)) denotes the nucleon state with four momentum p = (po,p) and spin s normal- 
ized covariantly as (N(p)\N(p')) = (2vr) 3 2p (5 3 (p - p'). U°^(q) gives the main contribution 
for n^ M (g) due to the perturbative contribution. On the other hand T^q) leads small 
contribution for n^ M (g), but the effect is vital contribution. 

If we consider sufficiently low nucleon density such as normal matter density (p^ ~ 0.17 
fm -3 ), the integral of the last term in Eq.( |2.2| ) can be approximated up to the first order 
of nucleon density p^ reasonably well. The linear density term corresponds to the matter 
with static nucleons (p = 0) and higher order correction terms correspond to the velocity- 
dependent terms involving the effect of Fermi motion (p ^ 0) and the complex interaction 
among nucleons. The linear expression can be calculated model-independently. On the other 
hand the higher order corrections depend on the model calculation, but in a few effective 
theory ]17| it is known that the effect for the linear result is fairly small (~ 10%) at nuclear 



matter saturation density. Hatsuda et al. also insist that the Fermi momentum correction is 
fairly small (~ 10%) up to twist-4 operators in fllQj . Thus we can safely neglect the effect at 
the saturation density. Therefore we can set p = (M N , 0) for T^ u (p, q), so that we proceed 
to discussions based on the assumption that all nucleons are at rest in nuclear matter. 

In Eq.( |2.2| ) the second term means the slight deviation from the properties in free state 
determined by 11°^. By applying QSR method to T^ u directly, we relate the scattering length 
extracted from the QSR for with the mass shift as one of deviation from the free state in 
the framework of QSR. Near the pole position of the J/ip, can be associated with the T 
matrix for the forward J/ip-N helicity amplitude ThH.h'H'^-, l \ where h(h') and H(H') are 
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the helicity of the initial (final) J/ip and the initial(final) nucleon, respectively. The relation 
between and Tkh^w is given by the relation 



^(9)^(0;, 9 )ef fc) (g) * ^ W^^). (2.4) 

l<? — -t- 2t j 



Here we introduce the coupling fj/^ and the J/^ mass mj/^ by the relation (0| J M | J/^ h \q)) 
— fj/tp m j/ip e ^ l K ( l) with the polarization vector normalized as J2h (l) e ^ (o) — ~ 
gV-v + q^q u /q 2 . Taking the spin average on both sides of Eq.( |2.4j ), T^ioj, q) is projected onto 
T(cu, q) = T^l (—3) and T hHh i H i(u, q) is projected onto the spin averaged J/ip-N T-matrix, 
T(u;, q). At low energy, q = (mj/^, 0) and p = (Mat, 0). T is reduced to the spin averaged 
J/ip-N scattering length aj/^ = l/3(2a 3 / 2 + ai/2) (04/2 an d 03/2 are the scattering lengths in 
the spin-1/2 and spin-3/2 channels, respectively) as T(mj/^, q = 0) = 8n(M N + rrij/^aj/^. 
We note that the negative aj/^p corresponds to attraction in our convention. 

We relate the parameters of the QCD Lagrangian with the hadronic mass and coupling 
using the dispersion relation. If one utilizes the retarded correlation function as a useful 
quantity for dispersion analysis, we obtain the following dispersion relation for T(u>, q); 

TV, 0) = I f *■ " {U,0) = - f * 2 4^4. (2-5) 

TT J -00 U — UJ — IE TT JO U — UJ 

Here the spectral function p(u, q = 0) is given with three unknown phenomenological pa- 
rameters a, b, c in terms of the spin-averaged J/ip-N forward T-matrix T such as 



p(u, q = 0) = — Im 

TT 



~/j/V> m J/V> 
( m2 - m j/i, + »e) 5 



T(u, 0) 



+ • ■ ■ (2.6) 



= a5'(u 2 - m]^) + b5(u 2 - m^) + c5(u 2 - s ). (2.7) 

• • • term in Eq.( |2.6| ) represents the continuum contribution and 5' in Eq.( |2.7p is the first 
derivative of 5 function with respect to u 2 . The first a-term is the double pole term cor- 
responding to the on-shell effect of T matrix and the coefficient is associated with the 
scattering length aj/^ as a = 8tt fj /^rrij /^{M N + mj /^)aj The second 6-term is the simple 
pole term corresponding to the off-shell effect of T matrix. The third c-term is the contin- 
uum term corresponding to other remaining effects, where Sq is regarded as the continuum 
threshold in vacuum. Now the contribution from the inelastic channels is not included in 
the ansatz of Eq.( [2.7|) . In this system the OZI rule restricts the inelastic channels of J/ip-N 
interactions to those containing charmed quarks, for example, J/ip + iV — > D + D + N and 
J/ijj + N — > A c + D. But all these processes are forbidden at the threshold. So fortunately 
this system is immune from such inelastic contributions. 

The parameters a, b and c in Eq. ( |2.7|) are not completely independent. That is, among 
these parameters we introduce a constraint relation, which is imposed by low energy theorem 
for the J/ip current-nucleon forward scattering amplitude. In the low energy limit uj — > 0, 
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T(u,0) become equivalent to Born term T Born (uj,0), which is zero in J/ip-N system. Now 
we get the following constraint relation from the low energy theorem, 



+ 



ni 



m 



+ - = 0. 

so 



:2.8i 



Therefore the spectral function is parametrized with two unknown phenomenological pa- 
rameters a and b by removing c from Eq. G2.8|) . The phenomenological (PH) side for II 
can be expressed as the combination of pole position for 11° and T^ u such as 



n 



NM 



9»v 



+ 



Pn 



+ 



m j/4> ~ (f 2M N [ {m 2 j^ - q 2 ) 2 m 2 ^ - q 



+ 



NM 
fiu 



OC 



F + AF 



m% + Am 2 , 



+ 



(2.9) 



j/i> 



where the pole residue F in IT^ is equivalent to fj/^Trij,^ and the deviation AF is pjqb/2M N . 
The quantity expressed as the shift of the squared J/ ij) mass in nuclear matter, 



Pn 



2M N f 2 , b m\ 



Pn 
2M, 



A? 



1%(M N + mj/^aj/^ (2.10) 



is proportional to the scattering length aj/^ through the double pole term in T^. Thus we 
can calculate the mass shift 5m jm in Eq. (|2.10| ) from aj/^ obtained by QSR for T^ v . 

We explicitly write down PH side with the unknown parameters a and b for T(q 2 ) using 
Eq. (|2.5| ), (|2.8| ) and ( |2.9| ). We take the n-th derivative with respect to q 2 after dividing T ph 
by q 2 as follows and define it as T^ 1 '. 



n\ \dq 2 



T ph (q 2 



f^ ph (g 2 ; a, b) 



m 



j/i> 



(n + l)m 2 j^ 
(m 2 ^ - q 2 ) n+1 



+ 



2\n+l 



(s - q 2 ) n+1 



1 



(«o - q 



(2.11) 



In order to construct the QSR we calculate the n-th derivative of OPE side similarly in 
the next section. 



3 The calculation of the Wilson coefficients for T uz/ 

Now we give the OPE expression for T^ v . The main task in OPE side is to calculate the 
Wilson coefficients based on perturbative QCD. In the case of J/ip the charmed quark mass 
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is so heavy that the calculation of the Wilson coefficients must be carried out explicitly with 
the effect of heavy quark mass. We now expand local operators up to dimension-4 in the 
OPE side. Then pure gluonic contributions must be only taken into account for the local 
operators. Up to this order of the OPE the nucleon matrix elements of two-gluon operators 
(GG) are most dominant. We note that contrary to vacuum QCD sum rule a new feature in 
T^y is that the nucleon matrix elements survive not only Lorentz scalar operators but also 
nonscalar operators. That is, we must consider new contributions from twist-2 operators 
with two spins for the Wilson coefficients. In order to calculate the coefficient function, we 
adopt the valid and well-known method for massive quarks propagating through the couple 
to soft gluons working as the external field, namely fixed-point gage method [18|. This gauge 
condition is expressed as x^A a „(x) = 0. The nucleon matrix element of two-gluon operators 
can be decomposed into the scalar part and the twist-2 part with an additional four-vector 
■ 2 - 1), through the simple tensor analysis ||19 . 



u„. ( u 



96 



(° 2 ) N (gorges - gasgp 7 ) -a((u-g) 2 - ^G 2 ^ 



- 2 (g aj upus - gasupUj - g^u a u & + gpsUaU^) 



— gasgp-/) 



(3.1) 



where we define (u ■ G) 2 = G a KX G a p A u K u p . By introducing the u^, one can imagine uniformly 
moving uncorrelated nucleons (p M = Mjyu M ) in nuclear matter, but in this case we set 
u = (1,0). The OPE expression for is written as follows by the combination of Eq.( |3.1| ) 
and the Wilson coefficients corresponding to each matrix element. 




f(n) OPE ^2) 



1 

3 



4 n) (0 



a: 



-G 2 

vr / n 



a 



—st(Gq u gi (Ji 

vr / n 



+ {^ n) (0 - D?\i) - D™(£)} < -ST 



)(") 



(»)/ 



'Oct, 



(3.2) 



Here we define the dimensionless parameter as £ = —q 2 /4m 2 (m c ; charmed quark mass) 
and p = £/(l + £)• ST means making the twist-2 operators symmetric and traceless in its 
Lorentz indices. In Eq.( |3.2| ) each coefficient function is given using Gauss hypergeometric 
functions 2-F1 for arbitrary q 2 as follows: 



4 n) (£) = - 2 " (W (2 t + 5);! (4^r ( " +2) (l + 0- ( " +2) 2*1 (n + 2,- 1 -, n+ 7 -;p) (3.3) 
D (n) {0 = 2" +3 (n + l)(n + l)! (4m , ) _ (n+2)(1 + _ (n+2) 



3(2ra + 3)!! 
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2(n + 2) „ / 1 5 



3(n + 2) 2 _ / 1 7 

+ (l + 0(2n + 5) »*H W + 3 '2' W+ 2 ; ' 

2^(n + l)(n + 2)! (n+2) (n+2) 

3(2n + 5)!! 1 J 1 +U 



(3.4) 



x 



2 Fi (n + 2, J, n + 



n + 2 / 1 7 

a Fi n + 3, -,n + -;p 



2(1 + 



(3.5) 



)(n)m _ 2"+ 3 (n + l)(n + l)! (i ^ ) _ (TO+2)(i | _ (n+2) 



x 



3(2n + 5)!! 

n + 2) 2 F 1 fn + ^.n + ^p) +4(2n + 5) 2 F X (n + 2,^,n + ^;p) 



(3.6) 



The Wilson coefficient of Eq. (p.3|) for scalar operator have already given in [5 i] and Eq. (|3.4j) , 
( |3 . 5| ) and (|3.6|) are new contributions for twist-2 operator. Eventually by equating Eq. ( p.ll] ) 
and Eq. (|3.2|) we obtain the moment sum rule expressed as the form of the n-th derivative 
with respect to q 2 , 



f (n)ph (e; a,b) =f (n)OPE (0- 



(3.7) 



The manipulation of the derivative ensures the enhancement of low energy part not to 
depend on the details of high energy part. The vacuum sum rules have been utilized for 
investigation into the free state of charmonium by Reinders et al. [Ill in moment sum rule 



and by Bertlmann [^UJ in Borel sum rule. Furnstahl et al. have studied the spectra 
of J /if) at finite temperature using both QCD sum rules. We summarize the well-known 
behavior of the moment sum rule for the variations of n and q 2 . 

• The convergence of OPE side is worse with n larger but better with q 2 larger. 

• In contradiction to this behavior of OPE side, the unwelcome contributions from the 
continuum in PH side grow with q 2 larger but decrease with n larger. 

We must choose the reliable stability region of moment sum rule for the change of the both 
n and q 2 . The moment sum rule for T^ v is the method to investigate the deviation from the 
properties in vacuum obtained from IT . So we should adopt the same regions of n and q 2 
as Ho, which can reproduce the nature of J /if) in the moment sum rule reasonably well. 
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4 Moment sum rule for n 



o 



In this section we calculate the window of n for various values of £ by applying the moment 
sum rule to vacuum correlation function Hq [|TT] . 



In OPE side ra-th derivative for IIq is expressed as 

\0 



where cj n) (0, cT(0 

are given by 



1 (J_ 

n\ \dq 2 

1 

3 



lL 



il 



in). 



a 



s ri2 



71 



(4.1) 



9 2 n (n + l)(n - 1)!, 
4^2 (2n + 3)!! 

(2n + 1)!! 2n + 3 



(4m 2 c )- n (l + 02 n Fi [n,-,n+-;p 



(4.2) 



3.2-%!2(n + l) 2 F 1 (n,| ) n+f;p 

1 













X 


7T — 






(I" 








V 2 





3 

47T 



n 



2.F1 (n, l,n + 2;p) 



+ 



1 



3(n + l)(n + 2) V2 4vr 



TT 



2-Fi (n,2, n + 3;p) 



vr_^_\_ 2 ln(2 + Q 2 + £ 2^(^+1,1,71+1^ 



47T 



7T 



(4.3) 



In Eq. (|4.1|) the OPE side of moment sum rule in vacuum n| ) n ' ) is the combination of the first 
term from a bare loop contribution and the second term replaced by the expectation values 
of the nucleon into it of the vacuum. On the other hand the relation between the J / ip mass 
of lowest resonance and IIo in PH side is given as 



ni n) (0 = 7 



9 m 



J/4> 



1 



4 9j /4) {mj^ - w 1 



(4.4) 



where mj/^ and gj/^ are the parameters of the lowest-lying resonance and S n represents 
the contributions from higher resonances. Moreover in Uq 1 ^ we add continuum contribution 
proportional to 1/Att x (1 + a s /ir), which we for simplicity assume to be constant value 
without the dependence of charmed quark mass. From Eq. (|4.4j) we eliminate the coupling 
parameter gj/^ by taking ratios of the n-th moments and the [n — l)-th moments. Finally 
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the bare mass of J /if) is derived from the following relation, 



m 



co 2 + 



nr i} - £ (i + f ) 


i i 


n-1 (sq-uj 2 )™- 1 




) 1 1 



1/2 



(4.5) 



First we fix £ to be from 0.0 to 3.0 at 0.5 intervals. These values are equivalent to the 
magnitude of from to 4 [GeV] in \J —q 2 . In Fig.l we show the J /if) bare mass determined 
for the change of £. Here we used sq = 3.6 2 [GeV 2 ] adopted in [[□]]. We must 



from Eq.( 

read off the range of n for each £ from the Figure 1. The window of n corresponds to find 
the region stabilizing the J / if> bare mass for the change of n. Then we obtain the windows of 
n for each £ as follows: n x = 2, 3, 4 for £ = 0.0 (-g 2 = 0.00 [GeV 2 ]), n 2 = 3, 4, 5 for £ = 0.5 
(-g 2 = 3.18 [GeV 2 ]), n 3 = 4, 5, 6, 7 for £ = 1.0 (-g 2 = 6.25 [GeV 2 ]), n 4 = 5, 6, 7, 8 for £ = 1.5 
(-q 2 = 9.23 [GeV 2 ]), vn 5 = 6,7,8,9 for £ = 2.0 (-q 2 = 12.1 [GeV 2 ]), n 4 = 7,8,9,10 for 
£ = 2.5 (-g 2 = 14.9 [GeV 2 ]), n 5 = 8,9, 10, 11 for £ = 3.0 (-q 2 = 17.6 [GeV 2 ]). These points 
seem to reproduce the bare mass reasonably well for the experimental value mj/^ = 3.096 
[GeV]. 



5 Numerical results 



By inserting the sets of £ and n obtained in section 4 into Eq. fl3.7| ), we can determine 
unknown parameters a and h simultaneously by fitting the left-hand side to the right-hand 
side. Concretely the order of calculation is as follows; At first we arbitrarily choose two 
points in the window of n for the fixed £ and make a simultaneous equation for a and b by 
inserting two n chosen. We consider all such combinations for each £ and take the average 
of a solved for each combination. Eventually the scattering length is easily obtained from 
a. 

To calculate the scattering length we use the next values for other various parameters. 
First of all, in PH side we adopt m j/ip = 3.1 [GeV], M N = 0.94 [GeV] and s = 3.6 2 [GeV 2 ] 
adopted in 11]. The coupling is determined from the experimental value of Tj/^ as follows: 



opee 

61 JH> 



j/i> 



47re 2 a 2 m J/v , 



1.7 x 10" 



(5.1) 



where e q is electric charge of quark (e c = 2/3 for charm quark) and a is the fine structure 
constant (= 1/137). Indeed this values of the coupling can be also determined by sub- 
stituting the experimental values of J/ip bare mass into Eq.( |4.5| ) oppositely. Then QSR 
for I! reproduces well the experimental values of the coupling [fll}l . For QCD Lagrangian 
parameters we use the following functions dependent on £ given in 

a s (4m 2 ) 
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1 + & a s {Aml) ln(l + 



a s (4m c 



0.3 



(5.2) 



m c (f) = 1.28 x 



ln(2 + 0-2 hi2 



GeV 



(5.3) 
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0.0 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


- ay alar [fm] 


0.091 


0.068 


0.070 


0.063 


0.059 


0.057 


0.055 


5m scalar [ Me y] 


5.3 


3.9 


4.0 


3.6 


3.4 


3.3 


3.2 


- a^ lst -' 2 [fm] 


0.120 


0.090 


0.092 


0.083 


0.078 


0.075 


0.073 


6m^ lst - 2 [MeV] 


6.9 


5.2 


5.3 


4.8 


4.5 


4.3 


4.2 



Table 1: J/ip-N scattering length and the mass shift of J/ip in the case of the scalar operator 
only and the scalar + twist-2 operator at normal matter density pn = 0.17 [fm -3 ]. 



In OPE side we determine the nucleon matrix elements as follows: 

'0£s 



G 

7T / JV 



7T 



Oct, 



N 



•(1.222 ±0.282) [GeV 2 



(0.094 ±0.010) [GeV 2 



(5.4) 



(5.5) 



The scalar part is evaluated from the ir-N sigma term [19, 22]. The twist-2 part is deter- 
mined from the gluon distribution function of a nucleon, which is obtained by leading order 
parametrization to the experimental data of deep inelastic scattering |19|, ^3]. 

We list the results in the case of scalar operator only and the case involving twist-2 
operator in Table 1. 



6 Concluding remarks 

The direct application of moment sum rule to the forward J/ip-N scattering amplitude 
supplies us the fascinating result for the J/ip-N interaction. That is, the J / tp-N scattering 
length a j indicates negative value (about —0.1 fm). This result suggests that the attractive 
J/ijj-N interaction is not sufficient to form a bound state with one nucleon, but it could 
make a bound state with nuclei. The absolute value is certainly smaller than the typical 
hadronic size 1 fm and the scattering length of light vector meson- N systems (a p ~ —0.47, 
~ —0.41, ~ —0.15) 0, but the experimental creation of J/ip at the threshold would 
lead to formation of a bound state inside a heavy nucleus. Our result is smaller than those 
obtained recently by S.J.Brodsky et al. pi] and G.F.de Teramond et al. in the QCD 



sum rule approach . Their method is based on the on-shell calculation on charmed quark 
mass (q 2 = 0). 

In this study we newly calculated the Wilson coefficients for twist-2 gluon operators 
(dimension-4) in the form with quark mass explicitly. The nucleon matrix elements of 
twist-2 gluon operators is about 1/10 times as large as that of the scalar part, but the total 
contribution with the Wilson coefficient makes the absolute value of scalar part larger about 
30%. 
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From aj/^, we can estimate a total cross section (aj/^ = Aircij,^). The result is about 
1.26 mb at the threshold. The contribution from the elastic channel corresponds to 20 % of 
the nuclear absorption cross section derived from the experimental data, a na b = 7.3 mb. 

Next in the linear density approximation we can calculate J/ip mass shift from aj/,^. 
The result gives very small decrease of mass (about —4 to —7 MeV), about 0.1 to 0.2% at 
normal matter density. Since the slight mass shift is of the order of MeV, the change is 
sufficiently larger than the leptonic decay width of the order keV. So we consider J/ip is a 
good probe for the observation of medium effect. 

Recently just before I submit this paper, Klingl et al. reported the mass shift of about 
—20 MeV (0.6%) for J/xjj in nuclear matter [27]. The pattern of derivation is similar to the 
case for light vector meson & 0. 
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Figure Captions 



Fig. 1 Stability region in n for various values of £. For comparison the experimental mass 
values have also been indicated. 

Fig. 2 J/ip-N scattering length determined for various values of £ and the stability region 
of n obtained from Fig.l. 
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